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We show that the convex set of separable mixed states of the 2x2 

■ system is a body of a constant height. This fact is used to prove that 
O || the probability to find a random state to be separable equals 2 times the 
! 1 1 ■ probability to find a random boundary state to be separable, provided 
£^ ' the random states are generated uniformly with respect to the Hilbert- 

Schmidt (Euclidean) measure. An analogous property holds for the set of 
positive-partial-transpose states for an arbitrary bipartite system. 

The phenomena of quantum entanglement became crucial in recent devel- 
. i-H , opment of quantum information processing. In general it is not easy to decide 

' whether a given mixed quantum state is entangled or separable pQ. The situa- 

tion gets simpler for the qubit-qubit and the qubit-qutrit systems, for which it 
is known that a state is separable if and only if its partial transpose is positive 
|21 EH- However, even in the simplest case of the two-qubit system, the struc- 
ture of the set of separable states is not fully understood. This 15-dimensional 
convex set M.g is known to contain the maximal ball inscribed in the set of 
all mixed states M.^ @]. Although some work has been done to estimate the 
volume of the set of separable states [5] El E| |H1 El and to describe its geometry 
^l^J, the exact volume of the set of separable states is unknown even in this 
simplest case [r^] . 

In order to elucidate properties of the set of separable states Slater studied 
the probability that a random state be separable, inside the set of mixed states 
and at its boundary, using the Hilbert-Schmidt measure. For the two-qubit 
system he found numerically that the ratio Q between these two probabilities 
is close to 2 ^3 |2J. The fact that it is more likely to find a separable state 
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inside the total set of mixed states than on its boundary is consistent with the 
existence of the separable ball centred at the maximally mixed state. 

We prove that for this system the ratio f2 is indeed equal to 2, provided we 
work in the Euclidean geometry and use the measure induced by the Hilbert- 
Schmidt distance. Furthermore, we show that the space of two-qubit (or qubit- 
qutrit) separable states can be decomposed into pyramids of the same height. 
These results hold also for the convex set of states with positive partial transpose 
(PPT states) for any K x M bipartite system. 

We will need some generalities about convex sets. Let X € MP be a convex 
body in W. D (in particular it is a closed, bounded set with nonempty interior) and 
dX its boundary. Let V denote the volume of X and A the (D — 1-dimensional) 
area of dX, both computed with respect to the standard Euclidean geometry. 
For each body we define the dimensionless ratio 

**> - ^ ■ o 

where r = r(X) is the radius of the maximal ball inscribed inside X. This ball 
is known as the insphere. Before we see why 'j(X) is interesting, we invoke 
some geometric definitions. A face of a convex body X is the intersection of 
X with an affine hyperplane which intersects dX, but does not separate (the 
interior of) X into two pieces. Such an affine plane is known as a supporting 
hyperplane. Such hyperplanes are important because any convex body equals 
the intersection of the half-spaces defined by its supporting hyperplanes. The 
polar Y° of a set Y C R D is defined as the set {z : (z, y) < 1 for all y e Y}. 
For additional background information consult the literature [L31HB] , 

We can now formulate a lemma which is stated for bodies whose insphere is 
the unit ball (i.e., the ball of unit radius centred at the origin). However, this 
hypothesis is just a matter of convenience and can always be achieved by proper 
re-centring and re-scaling. Note that while, in general, there may be more than 
one inscribed ball with maximal radius, this ambiguity cannot occur for bodies 
verifying the conditions of the Lemma. 

Lemma 1. For a convex body X whose insphere is the unit ball the following 
three conditions are equivalent: 

1. Every boundary point of X is contained in a face tangent to the inscribed 
ball. 

2. X = Y° , where Y is a closed subset of the unit sphere such that the convex 
hull of Y contains the origin in its interior. 

3. j(X) = A/V = D. 

Before we prove this, let us observe that the equality A/V — D, from condition 
3, is Archimedes' formula for the volume of a Z?-dimensional pyramid of unit 
height — where we define a pyramid as the convex hull of the base and the apex 
of some convex cone. Hence it will be true that 7(A) = D whenever the body 
can be decomposed into a union of disjoint pyramids of fixed height with apex 
at the center of the inscribed ball. This is true for regular simplices, cubes and, 
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more generally, for any polytope circumscribed around the unit ball, and hence 
verifying condition 1. The equality A/V = D holds also for balls, which can be 
explained by the ball being a limit of an infinite sequence of polytopes. Thus 
condition 1 can always be thought of as asserting such "decomposability into 
pyramids" , and we will refer to a body obeying it as a body of constant height. 
However, whenever the bases of pyramids (faces of X) are of dimension lower 
than D — 1, the limiting procedures needed to prove equivalences will require 
some care. (For the body of N x N density matrices D = N 2 — 1, while the 
maximal faces have dimension D — 2N +1.) 

Proof of Lemma 1: 1 =>■ 2: We first note that X equals the intersection of the 
half-spaces limited by supporting hyperplanes corresponding to faces tangent 
to the unit ball. Consider one such supporting hyperplanc. Since it contains a 
tangent face, it has a common point with the unit sphere. On the other hand, 
it cannot cut the unit ball in two, since it is supporting to X and the unit ball 
is contained in X. Hence the supporting hyperplane itself is tangent to the unit 
ball, and so the corresponding half-space is of the form {z : (z,y) < 1} for some 
y £ S ^ 1 . The set Y is the closure of the set of points y that arise this way. 
The condition on the convex hull of Y will be automatically satisfied as it is 
equivalent to boundedness of X. 

2 => 1: First, if Y is a subset of the unit sphere, then the unit ball is 
inscribed into the polar of Y with every y £ Y a contact point. A point lies 
on the boundary of X = Y° if and only if one of the inequalities that define X 
as a polar body becomes an equality, i.e., if it belongs to the face F y := {x £ 
X : (x, y) = 1} for some y £Y. F y is clearly contained in a hyperplane tangent 
to the unit ball and it contains the point of tangency (namely, y), so it is itself 
tangent to the unit ball. We have thus shown that every point of dX belongs 
to a face tangent to the unit ball, as required. 

1^3. With an insphere S ' 1 of radius 1 we can regard any one of its 
points cj £ S ^ 1 as a unit vector, and we can define a real valued function r(u)) 
such that r(u))co £ dX. Then 



where dco is the usual measure on the unit sphere. Similarly, if n(ui) is the unit 
normal of X at the point lying "above" u, 



Here, and below, we rely on the fact that the normal of a convex body is uniquely 
defined, and continuous, almost everywhere. Whenever the normal n(u)) is 
uniquely determined, the face containing r(uj)uj is also uniquely determined and 
hence perpendicular to n(u). By assumption 1, that face, and the supporting 
hyperplane containing it, are tangent to the unit ball. This leads to the equality 






(3) 



(r(oj)uj, n(u>)} = 1 



(4) 
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or {uj,n{uj)) = l/r(u>); substituting this into (3) and comparing with (2) we 
obtain V(X) = A(X)/D. 

3 => 1. We first note that (by an elementary argument using convexity) the 
set 3 := {uj e S ^ 1 : r(u)uj lies on a face of X tangent to the unit ball} is closed. 
If assumption 1 does not hold, the complement of 5 in S D_1 is nonempty and 
it is an open subset of this sphere, hence of positive measure. By the same 
argument as in the preceding proof we have, almost everywhere in S 15-1 \ H, 

(w,»(w))>l/r(w) (5) 
while still (ui,ti(uj)) = l/r(u>) for w£S. Hence, from (2) and (3), 

A{X) < [ duj r(u) D = DV(X) . (6) 

This concludes the proof. □ 

We can now easily establish 

Lemma 2. Any intersection of two bodies of constant height Xi,X2, containing 
the same inscribed sphere, is a body of constant height. 

Proof: Let Yj be the sets defined by X = Y° for j = 1,2. Then X x n X 2 = 
(Yi U Y2) and of course Y\ U Yz is a subset of the sphere, if each of the Yj's 
was. Making use of property 2 we conclude that X\ n X 2 is a body of constant 
height. □ 

We are going to use the above geometric concepts to investigate the set 
M.( N ) of density matrices acting in TV-dimensional Hilbert space. An operator 
p : TLn — ► T~In belongs to if it is Hcrmitian, p = p\ is (semi) positive 

definite, p > 0, and is normalized, Trp = 1. Due to the latter normalization 
condition the set M^ N > is compact and has dimensionality D = N 2 — 1. The 
sphere inscribed in MS N ^ is centred at the maximally mixed state = 1/N 
and has radius r = 1/ \J (N — 1)N, if computed with respect to the Hilbert- 
Schmidt distance, d 2 is (p, a) = Tr(p — a) 2 . Observe that r is equal to the radius 
of the sphere inscribed inside the (JV — l)-dimensional simplex An-i with edge 
lengths 

A recent analysis |E| of the volume and area of the set of mixed states shows 
that their ratio A/V is equal to y/N(N - 1) (N 2 - 1). Therefore the coefficient 
7 reads then 

A 1 



l{M {N) ) = r-= - — ^/N(N-1){N 2 -1)=N 2 -1 = D. (7) 

V y/(N-l)N 

This observation inspired us to propose 

Proposition 1. The set Ai 1 - 1 ^^ of mixed quantum states is a body of constant 
height. 

Indeed, it is immediate that condition 1 of Lemma 1 holds. A density matrix 
lies on the boundary of My** if and only if it has a zero eigenvalue. Such 
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a matrix belongs to the face consisting of all density matrices with support 
in the subspace of the Hilbcrt space that is orthogonal to the corresponding 
eigenvector. [All such faces are maximal and isometric to Ai^ -1 ^.] The point 
of tangency of that face to the inscribed (Hilbert-Schmidt) ball will be a density 
matrix with one eigenvalue (corresponding to the same eigenvector as above) 
zero and all other eigenvalues equal. [Like the set of pure states, both the set 
of maximal faces and the set Y of points of tangency have the structure of 
a complex projective space.] Of course, we can also derive the conclusion of 
Proposition 1 from J7J and the condition 3 of Lemma 1, but such an argument 
would not be self-contained and, moreover, would obscure the matter. 

Let us now take the dimension N as a composite number, say N = KM, 
and consider operators acting on a composite Hilbert space TLn = % (g> %. 
Such a decomposition of the Hilbert space allows one to define separable states, 
as those represented as a convex sum of product states |18| . 

L 

Psep = H Pf ® Pf ' ( 8 ) 

3=1 

where operators p A and p B act on Hilbcrt spaces Tin and TLm, respectively, 
the weights are positive, qj > 0, and sum to unity, Y2j=i 1j = 1. A state 
which is not separable is called entangled. Any separable state has a positive 
partial transpose, Ta(p) = (T <8> t)p > 0, and this criterion is sufficient if 
iV = 2x2 = 4or7V = 2x3 = 60|3]. Thus the sets of separable states and 
PPT states (states with positive partial transpose) coincide in these cases and 
are equal to the intersection of and its image Ta(M,( n >) - see Fig. ^ In 

any dimension N = KM the sets A4^ N ' and Ta(M^ n ') have the same shape, 
volume and surface area, because the partial transpose Ta acts as a reflection 
with respect to the affine subspace of PPT-invariant states (which includes the 
maximally mixed state p„). 




•tot 



Figure 1: The set M P p T of PPT states arises as a common part of the set of 
mixed states MS- N ' and its image Ta(MS n ') with respect to partial transpose. 
Its volume and area are denoted by Vppx and A-p-pT- 
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Let us denote the volume and the area of the set of mixed states by Vtot = 
VoliM^) and A tot = Vol(9.M (Ar) ). Analogously Vp PT and A PPT represent 
the volume and the area of the set of PPT states. The key result of this work 
will follow from 

Theorem 1. The set A4pQ T of PPT mixed quantum states is a body of 
constant height. 

Proof: The set M.( N ' of mixed states is a body of constant height, and so 
is its image, Ta (A4 ( N > ) . Both sets contain the same inscribed sphere of radius 
r = 1/ y (N — 1)N. Since the set of PPT states is an intersection of two bodies of 
constant height, M^pp T = fl Ta{M^ n " > ), the intersection lemma (Lemma 

2) implies that this set has constant height. □ 

The total area of the set Ai^ N ' of mixed states can be divided into two parts: 
one part Ap containing PPT states, and another part A^p containing not PPT 
states - see Fig. By definition A to t — Amp + Ap. On the other hand, the 
surface of the set of PPT states is a union of two congruent parts: the part Ap 
described above, and its isometric image under Ta - see Fig. \T]p. Accordingly, 
it is possible to infer that the area ^4ppt of the set A4p PT equals 2Ap once we 
establish that the area of the intersection of the two parts (the "corners" of the 
set PPT on Fig. ^) is zero. This is indeed the case and is the content of the 
following 

Lemma 3. The D — 1-dimensional area of the intersection of the boundary 
dM.^ with its image under Ta is zero. 

Before we prove this, we need to recall some basic facts about the boundary 
dM {N K First of all any density matrix can be written as 

p = UEU f , (9) 

where U is a unitary matrix and E is diagonal. As indicated in the comments 
following Proposition 1, dAi^ N ' consists exactly of those density matrices that 
have a zero eigenvalue. Further, after disregarding the density matrices with 
multiple eigenvalues (a subset of dimension D — 2, and hence of surface measure 
zero), the remaining "generic" part of N > is a smooth manifold which — 
because of eq. — is naturally diffeomorphic to the Cartesian product of the 
open simplex A := {(0, A2, . . . , Ajv) : < A2 < . . . < Ajv, 52 J= 2 A; = 1} an< ^ 
the flag manifold Fl N — U(N)/U(1) N . The first component is just the ordered 
sequence of eigenvalues of the matrix, while the second component is the set of 
of unitary matrices with those matrices that commute with E "divided out" ; the 
second component is in fact given by the corresponding sequence of eigenspaces. 
(A more detailed analysis of the stratification of can be found in 

Moreover, as is well-known (cf., e.g., |17p. this diffcomorphism transforms the 
surface measure to a product measure of the type /(A) dX <g> du, where dX is the 
Lebesgue measure on A, /(A) - a positive continuous function on A, and du is 
the invariant measure on Fl N induced by the action of U(N). 

The heuristic idea behind the proof that follows is that the area of intersec- 
tion is zero if the normals of the two hypersurfaces are distinct in "nearly all" 
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places where they intersect. However, some care is needed in order to ensure 
that this argument works. 

Proof of Lemma 3: We will first prove a weaker statement, namely that the 
intersection dM.^ N ^ H T^fdM.^) has empty interior (relative to dA4^). 

We shall argue by contradiction. Assume that the intersection does have 
nonempty interior. Then the subset C gen of the intersection consisting of density 
matrices p verifying 

• p and Ta(p) belong to the generic part of 

• the pure states corresponding to the eigenvectors of p and Ta(p) are not 
separable 

also has a nonempty interior since each of the above two conditions is satisfied 
outside a (closed) set of dimension < D — 1. Let po be an interior point of 
Cgen ■ We now recall that the first eigenvector (of eigenvalue zero) of a generic 
boundary density matrix p corresponds to the pure state which is normal to 
dM.( N ^ at p. Thus, by our assumption, the normal to dM.^ at Ta(po) corre- 
sponds to an entangled state and, consequently, its image under T4 - which is 
the normal to Ta(0A4^ ns> ) at po - corresponds to a trace one matrix which is not 
positive-definite. (This is because for all entangled pure states their negativity 
is positive 12131^, so the PPT criterion is necessary and sufficient for separa- 
bility of pure states). In particular, the (unique) normals to TA(dM.^ N ') and 
dM.( N ) at po are distinct. This implies that, in sufficiently small neighborhoods 
of po, the intersection of Ta^-M^) and dAi^ N ' is a D — 2-dimensional object 
and so po cannot be an interior point of the intersection and a fortiori of C gen , 
a contradiction. 

Since dAi^ n Ta(9A^^) is closed, the above argument does shows that 
it is a rather "thin" set ("nowhere dense" in topological parlance). However, 
this is not sufficient for our purposes because intersections of convex surfaces 
may be rather complicated, in particular they may be of positive area in spite of 
having empty (relative) interior. To correct this shortcoming, we need to choose 
Po G Cg en to be a Lebesgue point of (the indicator function of) C gen rather than 
an interior point of Cg en - Then, for every sufficiently regular sequence (Nk) of 
neighborhoods of po (relative to dA4^) which shrinks to po, the ratios of 
the areas of Nk fl C ge n and of Nk tend to 1 (this is just a reformulation of the 
definition of a Lebesgue point). By the Lebesgue differentiation theorem, almost 
all points of a measurable set have this property (see, e.g., Theorem 3.21). 
Accordingly, if we assume that dM m n T A (dM w ), hence C gen , has strictly 
positive area, then such a choice of po G C gen is possible and we are led to a 
contradiction in the same way as earlier. □ 

The ratio between the probabilities of finding a PPT state in the interior of 
and at its boundary is defined by |13|. 

^ _ PV Vppx/Vtot _ Vppt A iot 

~ p A '~ Ap/Atot ~ V tot A P 
The above quantity allows us to formulate the main result of the work 
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Theorem 2. For any bipartite K X M system the ratio is equal to 2. 

Proof: By Theorem 1 the set of PPT states has constant height, hence ^(M < ppr > ) 
rAppT /Vppt — D. For any fixed system of size N = KM, the sets and 
M ( p N p\ have the same dimensionality D and the same radius r of the inscribed 
ball, so that 

^tot _ ^4ppt 

Vtot VppT 

Lemma 3 implies Ap = Appx/2. Substituting this relation to (|10|l and making 
use of we find 

n= V ^ A ™ = 2 . (12) 

Vtot Appt/2 

which concludes the proof. □ 

The relation SI = 2 is thus established for any bipartite system. In the 
simplest cases of 2 x 2 and 2x3 systems any PPT state is separable 3^ Hence, 
for such systems, eq. (|12fl describes the ratio between the probabilities of finding 
a separable state inside the set of mixed states and at its boundary. 

The value = 2 is consistent with the numerical data obtained in H^j 
for the Hilbert-Schmidt measure. Since our reasoning hinges directly on the 
Euclidean geometry, it does not allow to predict any values of analogous ratios 
computed with respect to Bures measure |22|. nor other measures. On the 
other hand, the result obtained is valid also for the set of real density matrices, 
for which the complex flag manifold Fl% has to be replaced by the real flag 
manifold, Fl* = 0(N)/0(1) N . The set of real density matrices is often used 
as an attractive toy model, since its dimensionality Dr = N(N + l)/2 — 1 is 
smaller than D = N 2 — 1 of the full set of complex states. 

The geometry of the 15-dimensional set of separable states of two qubits is 
not easy to describe. In this work we have established a concrete rigorous result 
in this direction, proving that this set is "pyramid-decomposable" and hence is a 
body of constant height. This is also true for the set of positive-partial-transpose 
states for an arbitrary bipartite system. We hope that these properties will be 
useful for further investigation of the geometry of the set of separable states. 
Although in this work we have concentrated on the bipartite case only, one 
could try to obtain similar results for a general class of multipartite systems, 
for which some estimates of the volume of the set of separable states or PPT 
states are known El ESI EU 
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